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A model of arsenic clustering in silicon is proposed and analyzed. The main feature of the 
proposed model is the assumption that negatively charged arsenic complexes play a dominant role 
in the clustering process. To confirm this assumption, electron density and concentration of impurity 
atoms incorporated into the clusters are calculated as functions of the total arsenic concentration 
at a temperature of 1050 °C. A number of the negatively charged clusters incorporating a point 
defect and one or more arsenic atoms (DAsi) - , (DAsi) 2- , (DAs2)~, (DAs2) 2 ~, (DAS3) - , (DAS3) 2- , 
(DAS4) - , and (DAS4) 2- are investigated. It is shown that for the doubly negatively charged clusters 
or for clusters incorporating more than one arsenic atom the electron density reaches a maximum 
value and then monotonically and slowly decreases as total arsenic concentration increases. In 
the case of cluster (DAS2) 2- , the calculated electron density agrees well with the experimental 
data. Agreement with the experiment confirms the conclusion that two arsenic atoms participate 
in the cluster formation. Among all present models, the proposed model of clustering by formation 
of (DAs2) 2_ gives the best fit to the experimental data and can be used in simulation of high 
concentration arsenic diffusion. 

PACS numbers: 61.72.Ji, 61.72.Tt, 61.72.Ss, 85.40.Ry 
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I. INTRODUCTION 

Using low energy high dose arsenic ion implantation, 
one can produce the active regions of modern integrated 
microcircuits characterized by very shallow junctions and 
high dopant concentrations. Thermal annealing is ap- 
plied after implantation for the arsenic activation and 
damage reduction. During the initial stage of anneal- 
ing, the arsenic atoms become electrically active, occu- 
pying the substitutional positions. Then a fraction of 
impurity atoms incorporates into the clusters, thereby 
decreasing the layer conduct ivityi^'^'£&Z'& & 1 1 1 1 Due 
to clustering, the concentration of charge carriers is less 
than the total impurity concentration because the clus- 
tered dopant atoms do not serve as a source of free charge 
carriers. We can observe the clustering phenomenon in 
many cases of the formation of highly doped semicon- 
ductor structures; for example, in the region with high 



impurity concentration during thermal diffusion of ar- 
senic atomsi 12 ' 13 ! 14 Deactivation of the electrically active 
dopant atoms also indicates that clustering can occur 
during thermal treatment of the supersaturated arsenic 
layers created by ion implantation with subsequent ther- 
mal or laser annealing i 4 ' 7 i 15 i 16 i 17 i 18 i 19 Moreover, a reverse 
annealing or transient increase of the carrier density can 
be observed if a doped specimen, previously annealed 
at temperature T\, is further annealed at higher tem- 
perature Tiw& It is supposed that transient dissolution 
of a part of the clusters occurs during such a reverse 
annealing, 20 

At very high dopant concentrations exceeding the solid 
solubility of As in Si, a significant fraction of arsenic 
atoms can form precipitates ,3iii£ i 6 i 8 i 10 i 15 i 21 i 22 The typi- 
cal distributions of the total dopant concentration and 
electron density after thermal annealing of silicon sub- 
strates heavily implanted by As were presented by No- 
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bili et aU (Fig. 7) and by SolmilS (Fig. lc). As 
can be seen from Fig. lcj*°- the carrier concentration 
reaches its saturation value n e = 3.57 x 10 8 /xm -3 , 
whereas the total arsenic concentration increases to the 
solubility limit C so i = 3.27 x 10 9 /im~ 3 — and contin- 
ues to increase in the region adjoint to the surface. It 
is supposed that the cluster formation occurs at high 
dopant concentrations j^'2' 3 -'£' 9 - ' 10 i 12 ' 15 i 16 more precisely, at 
the dopant concentrations approximately ranging from 
n e to 10 n e . If dopant concentration exceeds the solubil- 
ity limit C so i, precipitation occurs i 4 ** - Arsenic clustering 
and precipitation have attracted evergrowing attention 
of the researchers, as they are greatly important for the 
silicon integrated circuits technology. Both clusters and 
precipitates do not serve as a source of charge carriers and 
the total arsenic concentration may be higher than the 
electron density by about two orders of magnitude. The 
clusters and precipitates are metastable and can disso- 
ciate under subsequent thermal treatments. This means 
that stringent requirements should be placed upon the 
accuracy of the clustering models used to calculate the 
carrier concentrations in heavily doped silicon layers. 



II. ANALYSIS OF CLUSTERING PHENOMENA 

The development of a model for As clustering necessi- 
tates a consistent analysis of the experimental data and 
theoretical substantiation of the processes proceeding in 
heavily doped silicon layers. 



A. Clustering during thermal diffusion and 
annealing of ion-implanted layers 

A great number of clustering models with different 
types of clusters have been developed and proposed since 
1973^ 1 3 1 12 1 15 1 16 1 23 1 25 1 26 The variety of the pr0 posed mod- 
els reflects discrepancies in the experimental results 2 ^ and 
a great uncertainty in the nature of the clusters discussed. 
For example, inactive cluster VAs2 was considered in 12 
to explain the difference between the total and electri- 
cally active dopant concentrations for arsenic thermally 
diffused from a constant source. For the best fit to the 
diffusion profiles of ion-implanted arsenic, it was sug- 
gested that the electrically inactive arsenic atoms were 
incorporated into the clusters with four arsenic atoms per 
cluster X According to Tsai et ali 2 - who have also investi- 
gated diffusion in the ion-implanted layers, clustering of 
arsenic may be expressed as 
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where As + is the substitutional^ dissolved arsenic atom 
participating in the cluster formation; e~ is the electron. 

As follows from reaction (JTJ) , the clustered As atoms 
are electrically active at annealing temperatures, but are 
neutral at room temperature. 



It was found in 8 - that the formation of clusters ( As2 V) ; 
occurs due to the following reaction, which 
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•V 



(As 2 V) > 



(2) 



is most liable to cause deactivation because the con- 
centration of inactive As species increases to the fourth 
power of the active As concentration^ 1 * 1 ^ Here V 2- is the 
doubly negatively charged vacancy participating in the 
clustering. 

Indeed, if a local equilibrium is assumed, it follows 
from the mass action law for reaction ^ that 



c- 



AC Vx X 2 C 2 



A = H C, 



(3) 



(4) 
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where C A and C are the concentrations of the impurity 
atoms incorporated into clusters and substitutional^ dis- 
solved arsenic atoms, respectively; H A is the constant of 
the local equilibrium for reaction @; C' Vx and C l Vx are 
the actual and the equilibrium concentrations of neutral 

vacancies, respectively; x = — is the electron density 

. Ui 

normalized to the concentration of intrinsic charge carri- 
ers in a semiconductor during thermal treatment rij. 

One would expect from that reaction ^ yields the 
fourth power dependence C A oc C 4 , since for high con- 
centrations C » rij it might be reasonable to assume 
n C 

\ = — « — . However, as seen from Fig. 7 (Ref. 4 ) and 

m in 

Fig. lc (Refi* -), the electron density n « n e = const 
in the region of concentrations associated with active 
clustering. Consequently, reaction ([2|) yields the second 
power dependence C A ^ C 2 , which contradicts the ex- 
perimental dataiii 1 ^ 

Because of the vacancy participating in cluster for- 
mation, silicon self-interstitials were in turn ejected or 
left behind during clustering^ This assumption had been 
made much earlier in the stud y 15 i of impurity deacti- 
vation within the laser-annealed layers. However, the 
results obtained in* show that the point defects induced 
during clustering and/or precipitation make no contri- 
bution to the enhanced transient diffusion of arsenic im- 
planted at low energy, while these defects did coalesce to 
form extended defects at the projected range. Analysis 
of the defect microstructure has revealed that the tran- 
sition between arsenic clustering and SiAs precipitation 
is not abrupt, pointing to possible coexistence of arsenic 
clusters and SiAs precipitates (although precipitates are 
not directly observed) £ 
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The model 2 - was generalized by Solmi and Nobili 5 



lAs++e - Annealing , ro _ 1)+ 25^C ^ 



(0) 



to take into account a saturation behavior of the carrier 
density with increasing dopant concentration. For dif- 
fusion at a temperature of 1050 °C, the best agreement 
between the modeled and experimental curves describ- 
ing diffusivity against total arsenic concentration was 
achieved for m = 4. And at 900 °C, the neutral clus- 
ter model 12 provided a better agreement with the ex- 
perimental data. The reaction ([6]) was also used in 9 - to 
describe deactivation in the arsenic-doped layers, which 
were formed by ion implantation and high temperature 
annealing of silicon-on-insulator specimens. However, 
Solmi et al.— have recently used the following reaction 



toAs + + V 



As„ 



I 



(7) 



to model the transient arsenic diffusion. Here m is 
assumed to have the values between 2 and 4 to take 
into account the fact that the As cluster was formed 
around a vacancy with the subsequent injection of the 
self-interstitial. 

It is important to note that reaction ((6|) was first the- 
oretically studied by Guerrero et al. in paper, 2 ^ which 
will be considered below. 



B. Clustering during deactivation of 
supersaturated arsenic layers formed by laser 
annealing 

Let us consider the clustering models used for the ex- 
planation of deactivation in the supersaturated arsenic 
layers formed by ion implantation and subsequent laser 
annealing. InAS it was found that a cluster consisting of 
four positively charged arsenic atoms and two additional 
negative charges gave the best least-squares fit with re- 
spect to the experimental data obtained in the case of 
deactivation within the supersaturated arsenic layers cre- 
ated by ion implantation and subsequent laser anneal- 
ing. If two negative charges are associated with a doubly 
charged vacancy, this model gains support from the the- 
oretical calculations of Pandey et ali^i Such a cluster has 
an intuitive geometrical and structural appeal — four 
arsenic atoms tetrahedrally arranged around a vacancy 
may relax inward, thus relieving lattice strains due to the 
As size effect^ Using this cluster species, a single-step 
clustering reaction was postulated in 16 
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The assumption of self-interstitial ejection during clus- 
tering was supported by the experimental data obtained 



by Parisini et alJ£ who have observed that extended de- 
fects at the projected range were extrinsic or interstitial 
in nature. The reaction 



(Asa) 
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(As 2 v)+ + r 



(9) 



was proposed in^ only for the initial step of the deac- 
tivation process in arsenic implanted silicon specimens, 
first laser-annealed and then thermally annealed. Here 
I - is the negatively charged self-interstitial. The results 
obtained by Rousseau et al. 17 confirm the conclusion that 
arsenic is indeed deactivated by vacancies, with a concur- 
rent injection of self-interstitials. The latter follows from 
the observed significant enhancement of the diffusion in 
the buried boron layer underneath the As structure sub- 
jected to deactivation. Investigations of the enhanced dif- 
fusion in buried boron layers were continued in^ where 
it was proposed that small arsenic clusters of various sizes 
were formed around a vacancy during deactivation with 
injection of the associated interstitial into the bulk. 

Considering that the majority of the models described 
in this section are based essentially on the experimental 
investigations of the defect subsystem, we first concern 
ourselves with these experimental results. 



C. Experimental investigations of defect subsystem 

The methods of transmission electron spectroscopy 
(TEM)) 4 ' 6 ' 15 ' 22 ' 27 extended x-ray-absorption fine- 
structure (EXAFS)^! 2 ^ 2 ^ Rutherford backscattering 
(RBS)r^2i positron lifetime measurements^ 2 ^^^ 
combined with the measurements of electron momentum 
distributions^i^L^iiM 1 . and other methods 2 ^ are 
commonly used to investigate silicon specimens doped 
by As. For example, TEM observations show that 
As-lased samples are completely free of any visible 
defects^ On the other hand, TEM analysis reveals the 
presence of extended defects after subsequent thermal 
treatment ! 15 ! 27 For example, small {113} interstitial 
loops begin to appear beyond the plateau region of 
the As profile after low-temperature treatment^ After 
high-temperature thermal treatment (>850°C), a low 
density of very small As-related precipitates (about 2 nm 
in diameter) is observed. The precipitated As fraction 
cannot be responsible for the total amount of inactive 

It was established that EXAFS measurements pro- 
vided particular details about the local atomic structure 
around the dopant and confirmed the fact that As atoms 
are substitutional after laser annealing. 28,29 Note that, 
according to the x-ray standing-wave spectroscopy mea- 
surements, the As atoms remained at the substitutional 
positions even after 85% of the electrical activity has been 
lost due to thermal annealing of the laser melted layers^ 2 - 
On the other hand, a large amount of nonsubstitutional 
arsenic was detected by Rutherford back-scattering after 
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diffusion from TiSi2 source^ These defects are inter- 
preted as precipitates, probably formed due to the stress 
induced by TiSi2 layer. Indeed, when thermal annealing 
is carried out after the complete removal of titanium sili- 
cide (i.e., without stress), one can observe the activation 
of all the arsenic atoms diffused in the silicon substrate^ 
Thus, such full activation indicates the significant influ- 
ence of stress on the formation of precipitates. 

According to EXAFS measurements^ subsequent 
thermal annealing (350 - 750 °C) of the As-lased layers 
leads to the formation of VAs m clusters including up to 
7 ± 4 As atoms around the vacancy. However, at high 
temperatures it was observed that the number of the 
first neighbors of Si has tended back to 4 atoms. This is 
probably due to the precipitates observed by TEM over 
this temperature range. Quantitative analysis of EXAFS 
data for the ultra-low energy implanted layers after rapid 
or spike annealing 30 has revealed a site for As that was 
different from the pure substitutional one, suggesting the 
presence of clusters of As atoms coupled to vacancies. All 
the observed phenomena, namely, a low value of the to- 
tal coordination number in the annealed samples, small 
As-Si bond length, and the presence of As- As coordina- 
tions, may be explained by the assumption that cluster- 
ing of As ions with vacancies takes place to form VAs m 
structures. 

Ion channeling and RBS were used ir* 3 ^ to complement 
the EXAFS and electron microscopy results. In has been 
found that the deactivated As atom is displaced with re- 
spect to the lattice sites and, moreover, the average dis- 
placement is constant over the temperature range 450 - 
900 °C, being equal to 0.23 ± 0.06 A. Also, ir>2i it was 
proposed that the formation of larger clusters V n As m 
(n w 5 and m « 10) occurred starting from VAs2 or 
VAs4. The formation of these larger clusters is in agree- 
ment with a decrease in the first nearest neighbors of 
As. At higher temperatures, an increase in the number 
of the first nearest neighbors of As to 4 was observed 
by EXAFS measurements^ This means that the clus- 
ter species may differ at intermediate (350 -750°C) and 
higher temperatures. 

Very interesting information about clustering 
was obtained during the experimental investiga- 
tions of the doped layers by positron annihilation 
spectroscop y 33 ' 34 ' 37 ' 39 ! 40 ' 41 and by positron annihilation 
spectroscopy combined with the Hall-effect /resistivity 
measurements^ 9 - The defect subsystem was generally 
investigated before and after thermal treatment. To 
support these measurements, the positron lifetimes and 
core electron momentum distributions were calculated 
for different vacancy-donor complexes. 41 Lawther et 
alM- have investigated the doped layers melted by an 
excimer laser to obtain the profile with constant As 
concentration and a sharp fall-off at a depth of 0.2 
fim. Arsenic deactivation was initiated by annealing at 
750 °C for 15 s or by conventional thermal treatment for 
2 h. It was found that VAs m complexes with the average 
values of m greater than 2 caused arsenic deactivation in 



heavily doped Si. Myler et al. 34 have studied the silicon 
layers (with arsenic concentration 4xl0 8 /mi 3 ) fully 
activated by laser melting and subsequently annealed 
for 15 s at 500 and 750 °C. The changes in the positron 
annihilation spectra after the thermal treatment at 
500 or 750 °C were also attributed to the formation 
of As m Si4_ m V ac complexes. In this experiment, the 
impossibility of determining the number of impurity 
atoms incorporated in the impurity-vacancy complex 
was established. Ranki at al^ 9 - have studied the samples 
implanted by As during MBE growth at 450 °C. The 
samples were annealed in N2, either by RTA at 900 °C 
for 10-170 s or furnace annealing at 800 - 900 °C for 
two minutes. Based on the positron annihilation and 
Hall-effect/resistivity experiments, Ranki et al^ 9 - have 
concluded that the dominant defect in As-grown samples 
was VAS3. The measurements demonstrated high 
concentrations (above 10 5 /im -3 ) of VAS3 complexes in 
heavily-doped silicon. Moreover, larger V-As complexes, 
probably V2AS5, may occur together with VAS3 at high 
As concentrations. A relative amount of V2 AS5 increases 
with annealing. This cluster is even dominant after 
annealing at 800 °C. The VAS3 and V2AS5 complexes 
become unstable at 800 °C and 900 °C, respectively, then 
their concentrations decrease. Cluster reconstruction 
may occur during cooling , starting from VAsi to VAS2 
and then to VAS3. The clusters VAS2 and VAS3 may be 
in turn transformed to the V2AS5 complexes. Moreover, 
some VAs4 and other vacancy complexes may be also 
formed, but with much lower concentrations. 

High energy electron irradiation is commonly used in 
positron experiments for the generation of nonequilib- 
rium vacancies, mobile even at room temperature and 
liable to interact with the dopant atoms i 37 ' 38 ' 40 ' 41 For 
example, positron experiments^ 7 - were performed both 
with the As-grown samples doped by arsenic in concen- 
trations of 10 7 and 10 8 /itm -3 and samples subjected to 
2 MeV electron irradiation at 300 K. It was established 
that heavily As-doped silicon contained VAS3 complex 
as a native defect. Before irradiation, the concentration 
of VAS3 was equal to ~ 10 7 /1m" 3 at a doping level of 
10 8 /jm -3 . After electron irradiation, one can observe 
the pairs (VAsi) formed by a vacancy and a single impu- 
rity atom. It was demonstrated in 3 ^ that the migration 
of these VAsi defects started at around 450 K, leading 
to formation of VAS2 defects. And, in turn, these de- 
fects were transformed to VAS3 defects at 700 K. The 
VAS3 defects were stable at 700 K, representing the dom- 
inant vacancy-impurity cluster in heavily doped n-type 
Si at this temperature. The formation of larger VAS2 
and VAS3 complexes was significantly dependent on the 
As concentration. This interpretation of the experimen- 
tal data was confirmed by the results obtained in4°- In 
these experiments, the electron irradiated samples were 
annealed isochronally (30 min) at 300 -1220 K. Irradia- 
tion of the heavily doped Si samples has produced mainly 
the vacancy-donor pairs (VAsi) with a small concentra- 
tion of divacancies. Considering that, after irradiation, 
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VAsi concentration was much higher than the initial con- 
centration of VAS3, no signals from VAS3 were observed 
for the As-irradiated As-doped sample. From the core- 
region electron momentum distribution measurements, it 
was found*^ that defects in As-irradiated samples may 
be identified as VAsi, in samples annealed at 600 K - 
as VAs2, and in those annealed at 775 K — as VAS3. 
It is well known that a drastic drop in the conductivity 
occurs when heavily doped Si is annealed at tempera- 
tures between 400 °C and 500 °C. This deactivation of 
the dopants is partially reversible by annealing at 800 - 
1000 °C. Taking into account these experimental data, 
Ranki et al4£ have suggested that the formation and an- 
nealing of VAS3 at 700 and 1100 K, respectively, are re- 
sponsible for the observed behavior of the conductivity. 
The experiments conducted ir*2& were completed by stud- 
ies in 41 associated with thermal treatment up to 1220 K. 
It has been found that dissociation of VAS3 began at 1100 
K, and at 1220 K these defects were annealed away. 

Thus, the analyzed data show that complexing of As 
atoms with vacancies occurs in the layers heavily doped 
by arsenic. Based on the experimental data, various 
V m As„ clusters are possible but VAS2 and VAS3 are the 
most likely. Taking clustering into account, one can ex- 
plain the phenomenon of compensation at high doping 
levels. Nevertheless, some aspects are not clearly under- 
stood. For example, VAS3 starts to dissociate at 1100 K. 
This annealing temperature agrees well with the obser- 
vation that the number of Si first neighbors tends back 
to 4 atoms22 at temperatures higher than 750 °C. At the 
same time, clustering occurs at temperatures higher than 
750 °C as welli 1 ^ This means that such a high tempera- 
ture clustering is hardly explained by VAS3 complexes. 



D. First principles studies 

Along with various experimental investigations, a num- 
ber of the theoretical ab initio calculations have been 
performed to explain deactivation of arsenic atoms in sil- 
icon. From the calculations of Pandey et al.^l it follows 
that VAS4 complex including a vacancy surrounded by 
four arsenic atoms is energetically favored over both sub- 
stitutional, isolated As in Si and substitutional Si-As4 
configurations. This cluster is electrically inactive, be- 
ing responsible for arsenic deactivation and structural 
changes in heavily doped silicon. Larger defect clusters 
(e.g., V2AS4) should also form during annealing, whereas 
VAS4 is only the first step in the clustering process. These 
investigations were continued hij^ 3 - where general V n As m 
complexes have been considered. The formation energies 
of the vacancy-donor complexes VAsi, VAS2, VAS3, and 
VAs 4 were found to be 2.47, 0.82, -0.53, and -2.39 eV, re- 
spectively. Moreover, it was proposed that the complex 
VAS2 was mobile, as were the VAsi pairs. As these pairs 
moved, they reacted with other defects to form larger, 
immobile complexes. When considering the mobile VAS2 
complexes, Ramamoorthy and Pantelides^ 3 - have tried to 



explain the coupled diffusion phenomenona and cluster- 
ing. It was supposed that VAS2 and VAS3 were the dom- 
inant complexes in the deactivated specimens near the 
enhanccd-diffusion thrcshok h 46 ' 47 A high rate of As dif- 
fusion observed during rapid thermal annealing 46 was due 
to VAS2 migration over a very short time period (6 - 60 s), 
when no extensive clustering could occur. During their 
motion, VAS2 complexes reacted with other defects and 
formed larger, immobile complexes, which decreased the 
As diffusivity. 

Nevertheless, further investigation is necessary because 
some theoretical results disagree with the experimental 
data. For example, i n 44 i 45 it was stated that the forma- 
tion of VAS3 or VAS4 defects was exothermic, but, accord- 
ing to the first-principles calculations,— buildup of VAsi 
or VAs2 structures was endothermic. This means that 
even at low arsenic concentrations all impurity atoms 
must form VAS3 or VAS4 clusters after long-time treat- 
ment. In point of fact, at low dopant concentrations, one 
can observe intensive diffusion by means of VAsi pairs 
rather than clustering. 

Berding et al.^ and Berding and Sher— have used 
the electronic quasichemical formalism to calculate a 
free energy of various clusters, including AsSi4, VAS4, 
VAssSii, VAs2Si2- Here As Si4 is the arsenic atom in the 
substitutional position. The neutral cluster composed 
of threefold-coordinated second-neighbor arsenic atoms 
DP (2) was proposed by Chadi et ali 4 ^ as an alternative 
to exothermic VAS4. The clusters DP (2), V2AS6, and Sis 
were also included in the calculations performed in4£ In 
contrast to^ Berding et al.^ Berding and Sher— take 
into account the ionized states and entropy of the clus- 
ter formation. The entropy is unfavourable for the for- 
mation of a large defect complex such as VAS4. There- 
fore, the complete free-energy calculation is needed to 
determine the role of VAS4 in deactivation. Based on 
the complete free-energy calculation, it was found that 
cluster VAS4 was neutral, in agreement with the previ- 
ous findings^ Also, the energy of this complex (-1.62 
eV) was in a rough agreement with the value given in. ?4 
The VAS3 and VAS2 clusters were found to have one and 
two acceptor levels within the energy gap. Consequently, 
their formation energy may be effectively decreased when 
the Fermi energy is near the conduction-band edge. The 
equilibrium cluster concentrations depending on the tem- 
perature and concentration of dopant atoms are obtained 
using the minimized free energy of the system. For all 
the dopant concentrations and temperatures considered 
(400 - 1000 °C), three classes of clusters are dominant un- 
der equilibrium conditions: Sis, AsSi4, and VAS4. In all 
cases, the Fermi energy was near the conduction-band 
edge. Predominantly, VAS3 and VAS2 were singly and 
doubly ionized, respectively. At low arsenic concentra- 
tions (up to ~ 5 x 10 6 /um~ 3 ), noticeable deactivation 
was absent for temperatures higher than 500 °C. With 
the arsenic concentration raised to 5 x 10 7 /im~ 3 and 
higher, significant deactivation was predicted, mainly due 
to the formation of VAS4 clusters. Both concentration 
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and temperature influenced the contributions of different 
defects. For example, the concentration of VAS4 clusters 
has reached the concentration of isolated arsenic atoms 
in the lattice As Si4 for ~ 5 x 10 8 /im -3 at a temperature 
of 700 °C and for ~ 2 x 10 9 /Ltm~ 3 at a temperature of 
1000 °C. 

I n 25 i 48 the authors have attempted to explain the effect 
of the electron concentration saturation at high doping 
levels. However, the calculations performed in these pa- 
pers show that saturation is not reachable even at very 
high arsenic concentrations. For example, saturation was 
not reached with the arsenic concentration of 2 x 10 9 
/im~ 3 at a temperature of 700 °C, which is in disagree- 
ment with the experimental data (see Fig. 2from^.). It is 
important that, according to, 48 DP(2) clusters proposed 
by Chadi ct al. in 44 may be present as a deactivating 
species, insignificantly contributing to the deactivation 
under conditions of full equilibration. 

These new defects, first mentioned ini 4 - and called 
DP, represent a pair of two three-fold coordinated donor 
atoms. The lowest energy DP structures are DP(2) and 
DP (4), where the donor atoms occupy either second- or 
fourth- neighbor Si sites along the < 110 > direction. It 
was shown that DP (2) exist in the stable electrically ac- 
tive 2+ charge state that is donating electrons to the con- 
duction band or in the metastable neutral charge state 
capturing two electrons from the Fermi sea. A very im- 
portant conclusion made in^i is a threefold coordination 
of each dopant atom in DP(2) in the neutral trap state. 
Thus, relaxation of the neighboring Si atoms creates an 
increased volume, which is consistent with the experi- 
mental observations. Therefore, such an increased vol- 
ume around the donor atoms is possible without vacan- 
cies. The main ideas of 4 ^ were further developed in4^ 
Citrin et al4£ have proposed a new class of defects, called 
the donor-pair-vacancy-interstitials and composed of two 
dopant donor atoms near the displaced Si atom that is 
forming a vacancy-interstitial pair. This defect, which 
is like a hybrid of the donor-pair and Frenkel-pair de- 
fects, is denoted DP(i)V-I. For the case of Sb, the data 
of annular dark- field ( "Z-contrast" ) scanning transmis- 
sion electron microscopy (ADF-STEM) clearly demon- 
strate that, in heavily Sb-doped silicon grown at low 
temperatures, the primary deactivating defect contains 
only two Sb atoms*- 49 Thus, energetically favorable VSb3 
and VSb4 may be ruled out from the deactivation pro- 
cess. Only VSb2, DP (2), and DP (4) may be considered 
as liable candidates. Using additional STEM measure- 
ments, x-ray absorption data, and first-principles calcu- 
lations, Voyles et al42 have shown that neither VSb2 nor 
DP (2) and DP (4) defects are important in heavily Sb- 
doped Si. At that time, DP(2)V-I and DP(4)V-I were 
in line with the available experimental results, including 
positron annihilation spectroscopy data. 45 Both DP(2)V- 
I and DP(4)V-I are independent out of the pre-existing 
vacancy population, which is in conformity with the ob- 
servations of.— 

Nevertheless, as follows from the experimental data^ 



a deactivation mechanism can be different in the As- 
and Sb-doped silicon. The mechanism of arsenic clus- 
tering calls for further investigation. In,2S. Mueller et 
al. have investigated the electronic structure and charge 
states of various vacancy-impurity clusters using the 
first-principles density- functional theory (DFT). It was 
found that the VAsi complex can trap up to two con- 
duction electrons at a high n-doping level. Both VAs2 
and VAS3 act as a single-electron trap center. In agree- 
ment with,^ 5 -^^ Mueller et alM have found that in 
VAS4 the fifth valence electrons of all four As atoms are 
strongly bound and the complex is expected to remain 
neutral for any position of the Fermi level. Note that 
according to the calculations ofp 5 . the cluster VAS2 has 
two acceptor levels, whereas Solmi et al<22i have found 
from the carrier mobility measurements that the complex 
to be electrically neutral at room temperature. Ir* 5 ^ the 
first-principles calculations were performed to explain the 
coupled diffusion phcnomenona and clustering in heavily 
arsenic-doped silicon. As contrast to^ 3 - it was found that 
the VAS2 cluster is less mobile due to its high migration 
barrier. However, these clusters can contribute to the 
diffusion of As at elevated temperatures. 

Proceeding from the presented analysis, our under- 
standing of the clustering mechanisms is neither com- 
plete nor firmly established. This concerns the number 
of arsenic atoms in the cluster as well as its structure. 
However, a study based on the first-principles shows that 
clusters of a certain kind have an acceptor level and can 
be negatively charged. Let us consider the models for the 
clusters with different charge states. 

E. Models based on the mass action law 

The models based on thermodynamic formalism are 
usually used for simulation of high concentration dopant 
diffusion taking place during semiconductor processing. 
Assuming a local thermodynamic equilibrium among 
substitutionally dissolved arsenic, the dopant atoms in- 
corporated into clusters, and electrons, one can use the 
mass action law to calculate the concentration of clus- 
tered dopant atomsi2 i 12 ' 13 i 23 The charge conservation law 
for the reaction of cluster formation under the assump- 
tion of local charge neutrality is also useful for describing 
the clustering phenomenon. A widespread model of Tsai 
et al£ was the first to account for the saturation of the 
charge carriers at increasing doping levels.— It is interest- 
ing to compare this model with the latest experimental 
data. The mass action law for the reaction ([T]) at the 
diffusion temperature yields^ 

C cl =A cl nC 3 , (10) 

where C and n are the concentrations of clusters and 
electrons at the annealing temperature, respectively; A cl 
is the constant of local thermodynamic equilibrium. Af- 
ter cooling, all clusters become neutral and C ~ ur if 



7 



C ^> Uj. Here tir is the concentration of electrons at 
room temperature. It follows from the assumption of lo- 
cal charge neutrality that the electron concentration at 
the annealing temperature for C 3> rii is 



n w C + 2C* C ' = C + 2^ c 'nC 3 . 
Solving this equation, Tsai et ali^ obtained 

C 



1 - 2A cl C 3 



and 



C 1 =C 



3A cl C 4 



C + A cl C 4 



l-2A cl C 3 l-2A cl C s ' 
From Eq. (|13p it is inferred that under condition 



C — C sa t 



1 



(11) 



(12) 



(13) 



(14) 



both the concentration of arsenic atoms incorporated into 
clusters and the total concentration of arsenic atoms ap- 
proach infinity. The value nR e — C sa t is interpreted in 2 - 
as a maximum level of electron concentration. From the 
experimental data it is seen that a maximum equilibrium 
carrier concentration at the annealing temperature in the 
arsenic implanted layers may be described by the follow- 
ing relation^ 



1.3 x 10 n exp 



0A2eV 
k R T 



(15) 



Using Eqs. (JT3J) and (fig)), the value of A cl may be 
derived to be 



A cl = 



2n2" 



(16) 



at a temperature 



To illustrate, n e = 3.57 x 10 8 jtzm" 
of 1050 °C and A cl = 1.10 x 10~ 26 //m 9 . 

The values of ur as a function of the total arsenic con- 
centration calculated by the model of Tsai et ali£ are 
shown in Fig. 1. The experimental values of tir based 
on the data frorr£ is also given for comparison. As seen 
from Fig. 1, saturation of the electron concentration 
for the model 2 is observed at the extremely high total 
As concentrations, whereas the experimental data show 
that the carrier concentration reaches its maximum at 
C T w 8 x 10 8 /im" 3 . Thus, the model 2 is inconsistent 
with the experimental data. Moreover, as follows from 
the experimental data of Solmi and Nobili^ the value of 
n e obtained at room temperature corresponds to the pe- 
culiarity of the curve for diffusivity as a function of the to- 
tal arsenic concentration. Because this feature character- 
izes impurity diffusion at the annealing temperature, it 



was concluded in^ that n e had a physical meaning. Thus, 
the model 2 that includes cluster neutralization during 
cooling is also inconsistent with the experimental data 
for arsenic diffusivity. Indeed, the electron concentration 
at the annealing temperature calculated by this model 
and associated with the peculiarity of the diffusivity is 
three times greater than n e . 



10 a 



- model of Tsai et al. 

- neutral (VAs 2 ) 
■ neutral (VAsJ 

saturated value 




10° 



10 10 10 

Total arsenic concentration (nm~ 3 ) 

FIG. 1: Calculated electron concentration against the total 
arsenic concentration for different models of clustering: solid 
line — model of Tsai et al. (Ref .— ) , dash-dotted line — neutral 
clusters VAs2, and dashed line — neutral clusters VAS4. The 
experimental data (circles) are taken from Solmi and Nobil: 5 
for the diffusion at a temperature of 1050 °C. 

Fig. 1 presents the functions n = n(C T ) calculated 
for the cases when neutral clusters VAs2 and VAsa are 
formed. In these cases, saturation of the electron density 
is not observed and the calculations do not agree with 
the experimental data. The absence of saturation for the 
electron concentration corresponds to the first-principles 
calculations in^ (Fig. 2) and^ 8 - (Fig. 3). As seen from 
these figures, no saturation is observed at arsenic concen- 
trations when V Asi clustering is prevailing. 

The reactions 



mAs 



fee 



As! 



(17) 



were used in the basic model for the clustering of As in 
silicon proposed by Guerrero et al.— Here m and k are 
the numbers of arsenic atoms and electrons participating 
in the clustering process; r is the electric charge of a 
cluster. The function C = C(C T ) was investigated by 
means of the mass action law combined with the charge 
conservation law for the cluster formation reaction on the 
assumption of local charge neutrality 



= C + rC cl = C + rA ul C 



a. 



(18) 



It was supposed ir*22 that the charge of clusters be- 
comes zero when the sample is cooled down to room 
temperature. Therefore, C is interpreted as the concen- 
tration of electrically active arsenic at room temperature, 
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which equals the electron density n in the sample at room 
temperature. Four main conclusions have been drawn 
from the basic model: 

1. When the electrons do not participate in the re- 
action (fT7|) , i.e., k = 0, or the clusters are electrically 
neutral at high temperature, i.e., r = 0, C is a monoton- 
ically increasing function of C T and also there is no finite 
limit to C and electron density n at room temperature. 

2. If the clusters are positively charged at a high tem- 
perature (r > 1) and exactly one electron takes part in 
the reaction (JTSJ), i.e., k = 1, C is a monotonically in- 
creasing function of C T , which approaches a finite "sat- 
uration" value C max with increasing C T . 

As can be see from the reaction ([1]), the model of Tsai 
et al^ satisfies these conditions. On the other hand, from 
our calculations for this model (sec Fig. 1), it follows that 
saturation of the electron concentration is observed at ex- 
tremely high total As concentrations. It means that the 
saturation value obtained in such a manner is inconsis- 
tent with the experimental data. 

3. Provided that the clusters are positively charged 
at high temperatures (r > 1) and there are two or more 
electrons taking part in the reaction (fT7j) . i.e., k > 2, 
C increases with increasing C T only up to a finite value 
Cmax- A maximum concentration value C = C max is ob- 
served for C T — C max . With further increase in the total 
arsenic concentration C T > C max the concentration of 
the substitutional^ dissolved arsenic atoms C and hence 
electron density n at room temperature, are decreased. 

This behavior of n is at variance with the experimen- 
tal data, because experimentally measured n is approxi- 
mately constant with increasing C T (see Fig. 1). 

4. If the clusters are negatively charged at high tem- 
peratures (r < —1), then C = C(C T ) behaves as in case 
1, i.e., C, and hence electron density n at room tempera- 
ture, monotonically increase without limit for increasing 
C T . This behavior of n also conflicts with the experi- 
mental data. 

In what follows, the basic model is extended to take 
into account the different charge states of vacancies par- 
ticipating in the cluster formation^ It is shown that the 
main conclusions from the basic model are valid for the 
case when vacancy is involved in the formation of clus- 
ters. Thus, the generalized model of arsenic clustering 
proposed in2£ is not in agreement with the experimental 
data for all possible cases of the cluster formation con- 
sidered in this model. In our opinion, this contradiction 
is due to the assumption that charged clusters become 
neutral during cooling. As follows from^ the hypothe- 
sis that deactivation is taking place during cooling of the 
substrates contradicts the experimental data. 

Evidently, the models of arsenic clustering should give 
an adequate explanation not only for saturation of the 
electron density, but also for the plateau formed on the 
carrier concentration profile after annealing of the ion 
implanted layersjSJ^ As follows from the above analysis, 
there is no model satisfying these requirements. 

The performed analysis enables us to formulate the 



purpose of this study: the development of a more ade- 
quate model to explain the plateau phenomenon of the 
charge carrier profile at high concentration arsenic diffu- 
sion. 



III. MODEL 

To develop a model for arsenic clustering that can de- 
scribe the saturation phenomenon of the electron den- 
sity, we assume that not only neutral, but also negatively 
charged, arsenic clusters are formed. The possibility for 
the formation of the negatively charged arsenic clusters is 
confirmed by the calculations im 25 ' 26 It has been shown 
by Berding et ali2£ that clusters VAs^ and VAs2 have 
one and two acceptor levels, respectively. The calcula- 
tions carried out by Mueller et ali2£ also show that for 
the electron densities n > 10 4 fim~ 3 the majority of com- 
plexes VAs2 will be negatively charged. According toj2£ 
the complex VAS3 exhibits an acceptor level. This clus- 
ter is singly ionized at elevated Fermi levels, resulting in 
a loss of four mobile carriers for the formation of a VAS3 
aggregate. 

It is interesting that saturation of the electron den- 
sity is observed in the layers of GaAs heavily doped by 
silicon. To explain this phenomenon, it was supposed 
(see, for exampl o 51 i 52 ) that silicon atoms in heavily doped 
GaAs are dissolved substitutional^ not only in the Ga 
sublattice (donor SigJ, but also in the As sublattice (ac- 
ceptor Si^ s ). It was supposed in^ that a maximum elec- 
trical conductivity in Sb-doped silicon is limited by the 
defect formation at high dopant concentrations. The Sb- 
vacancy pairs were considered as defect dominating at 
high doping levels. Also, it was supposed that this de- 
fect acts as an electron acceptor and is responsible for 
saturation of charge carriers at high doping levels. Al- 
though a mechanism of such a saturation of the electron 
density in silicon heavily doped by Sb is different from 
that in As-doped silicon^ the idea of the limited elec- 
tron density due to compensation by the acceptors is a 
very fruitful one. For example, it has been supposed in — 
that phosphorus clusters are negatively charged. Due to 
this assumption, saturation of the electron density during 
high concentration phosphorus diffusion was explained. 

As distinct from the models of arsenic clustering con- 
sidered inj£i22 we assume that arsenic clusters have the 
same negative charge both at high diffusion temperatures 
and at room temperature. As follows from the analysis 
of the experimental data and theoretical investigations, 
the structure of arsenic clusters under diffusion temper- 
ature is still not conclusively established. Therefore, it 
is assumed that a point defect Di can be involved in the 
cluster formation and another defect D2 can be gener- 
ated during clustering. In this reaction for the 
formation and dissolution of clusters can be written as 



mAs + + miD ri + fce~< — >$ rci + m 2 D r2 , (19) 
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where A is the substitutionally dissolved arsenic atom 



As^ 



is the electron; $ is the cluster formed; mi and 



ni2 are respectively the numbers of defects Di and T) 2 
participating in the cluster formation; r\, r 2l and rci 
are the charge states of defect Di, defect D2, and cluster 
Q rci , respectively. 

The charge conservation law for the chemical reaction 
(fTi?)) is of the form 



m\Z\ - k = z C i+ m 2 z 2 , 



(20) 



where z\ and z 2 are the charges of defects Di and D2, 
respectively; zci is the charge of cluster <& rci in terms of 
the elementary charge. 

The mass action law for reaction (fT9|) yields 



rr 2 vm 2 



(21) 



where Kl is the rate of chemical reaction (|19|) in the 
forward direction; Kr is the rate of this reaction in the 
opposite direction; C is the concentration of impurity 
atoms participating in reaction (TT9"]) ; C r D \ is the concen- 
tration of defects in the charge state r\ participating in 
the cluster formation; C D 2 2 is the concentration of defects 
in the charge state r 2 generated during clustering; C$ is 
the concentration of the clusters in the charge state rci ; 
X is the electron density normalized to the concentration 
of the intrinsic charge carriers in a semiconductor during 
diffusion rii. 

Due to a high dopant concentration, the approximation 
of local charge neutrality can be used for the evaluation 
of x 



1 

2rii 



C + zci C<t 



+ J (C + zci Ct 



- 4nf 



(22) 



where Nb is the summarized concentration of acceptors; 
n ie is the effective concentration of the intrinsic charge 
carriers at a diffusion temperature calculated for a high 
doping level. 

It follows from equation (TT9")) that the concentration of 
impurity atoms incorporated into clusters C A = mC®. 
Moreover, considering highly mobile electrons, the mass 
action law for the reaction of defect conversion between 
the charge states is valid 



h T1 



-Zl 



C 



Dl 



Dl • 



C" 1 — h 

u m — n 



D2 



zi r<x — u 

°D2 — " 



-Z2 



DI ■ 



-z 2(J 



D2 



(23) 



(24) 



where and h r ^ 2 are the constants for the local ther- 
modynamic equilibrium in the reactions when neutral de- 
fects and are converted into charge states n and 
r 2 , respectively. 



Substituting (J23]) and ([24j) into (J2TJ) and taking into 
consideration that Cac — mCcii we obtain 



C AC =KC dX 



( k + r 2 m 2 - r-imx) qvh 



(25) 



where 



K 



K R (h D2r2 ) m 2 (C e «)™* 



mi 



(26) 



C 



D 



(C. 



Dl 



(Cd2 



\m 2 



Cdi 



f~<eq 
^DIX 



c 



D2 



C-D2X 
(~ieq 
D2 X 



(27) 

Here Cac is the concentration of clustered arsenic atoms 
and Ct — C + Cac is the total arsenic concentration. 
The parameter K has a constant value depending on the 
temperature of diffusion. This value can be extracted 
from the best fit to the experimental data. The quantities 



eq 

D1X 



and C D Q 2X are the equilibrium concentrations of 



C 

corresponding defects in the neutral charge state. 



IV. RESULTS OF CALCULATIONS 



Analysis of Eqs. (f2"2")l and (f2l)|) shows that saturation of 
the electron density will be observed for the formation of 
singly negatively charged clusters incorporating one ar- 
senic atom and a lattice defect. The calculations of the 
functions C = C(C T ), C A = C A (C T ), and n = n{C T ) 
for this case are presented in Fig. 2. For comparison with 
the experimental data^ the processing temperature was 
chosen to be 1050 °C (n» = 1.1665xl0 7 ^m~ 3 ). To cal- 
culate the functions C = C(C T ), C A (C T ), and n(C T ), 
a system of nonlinear equations (|22]) . (|25|) was solved 
numerically by Newton's method. In the case of singly 
negatively charged clusters (DAs)~ incorporating a lat- 
tice defect and one arsenic atom we adopt rci = -1 and 
77ii = 1- Eq. (|20[) for to = 1 yields that — zjtoi — z 2 m 2 + k 
= 2, and the system (|2"51). ([2^1 is of the form: 



c- 



KC DX 2 C 



(28) 



277,; 



C - KC DX 2 C - N E 



C-KC D x 2 C-N B ) +4n 



(29) 



As can be seen from Fig. 1, the concentration of charge 
carriers in the saturation region is approximately equal 
to the value n e = 3. 5665x10 s /7m~ 3 . In Fig. 2 for C T 
+00, the electron concentration reaches saturation at the 
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value n ~ 3.765xl0 8 /im -3 associated with the plateau 
of the charge carrier profile given in Fig. 1. The value of 
parameter KCd for this case was chosen as 1.0698 x 10 -4 
a. u. to provide a fit to the experimental value of n e . As 
seen from Fig. 2, the function n = n(C T ) is similar to 
that predicted by the model of Tsai et al.— Saturation of 
the electron density occurs at extremely high values of 
the total impurity concentration C T , which is in conflict 
with the experimental data. 
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FIG. 2: Calculated concentrations of substitutionally dis- 
solved arsenic atoms (solid line), clustered arsenic atoms 
(dashed line), and electron density (dotted line) against the 
total dopant concentration for the formation of singly nega- 
tively charged clusters incorporating one arsenic atom. The 
measured electron density (circles) is taken from Solmi and 
Nobili^ for the diffusion at a temperature of 1050 °C. 

In Fig. 3a the system of nonlinear equations (122p . (123)) 
is solved numerically for the formation of singly nega- 
tively charged clusters (DAs 2 )~ incorporating a lattice 
defect and two arsenic atoms (tqi = -1, mi = 1, and m 
= 2). Then, from (HOD, and {25} it follows that: 



C A =KC D x 3 C 2 



(30) 



X 



2n,i 



C--KC DX i C z ~N B 



C - -KC DX 3 C 2 - N B 



(31) 



The calculations for n = n(C T ) show that, with the 
increasing total concentration of impurity atoms C T , the 
concentration of charge carriers n reaches its maximum 
value n max =3.750xl0 8 /im" 3 at C^ ax - 1.634xl0 9 
/im~ 3 and then decreases monotonically. At maximum 
the concentration of the substitutionally dissolved ar- 
senic atoms C max is 7.943xl0 8 /urn" 3 . The parameter 



KCd used in this calculation is equal to 4.0xl0~ 14 /Ltm 3 . 
The calculated curve agrees with the experimental func- 
tion n = n(C T ) much better than those obtained by the 
model of Tsai et ali^ and the models of neutral (DAS2) 
and neutral (DAS4) clusters. 
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FIG. 3: Calculated concentrations of substitutionally dis- 
solved arsenic atoms (solid line), clustered arsenic atoms 
(dashed line), and electron density (dotted line) against the 
total dopant concentration for the formation of singly (a) and 
doubly (b) negatively charged clusters incorporating two ar- 
senic atoms. The measured electron density (circles) is taken 
from Solmi and Nobili^ for the diffusion at a temperature of 
1050 °C. 

The calculations of n = n(C T ) and C A = C A (C T ) 
were also performed for the clusters (DAs 2 ) 2 ~, (DAS3) - , 
(DAs 3 ) 2 ~, (DAs 4 )~, and (DAs 4 ) 2- . And similar behav- 
ior of n = n(C T ) was observed for all these cases. Also, 
the concentration of charge carriers increases with in- 
creasing C T , reaches a maximum, and then decreases 
monotonically. This decrease is more pronounced with 
an increased number of arsenic atoms in the cluster. Be- 
sides, the position of maximum carrier concentration is 
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shifted to smaller dopant concentrations with increases 
in the negative charge of the cluster. 

It follows from comparison of the calculated distribu- 
tions n{C T ) with the experimental data that the best 
fit is observed for the doubly negatively charged cluster 
(DAs 2 ) 2 " (see Fig. 3b). To obtain the curves shown in 
Fig. 4, we solve the system (|22p . (1231) having the follow- 
ing form for the clusters (DAs 2 ) 2 ~ : 



K C D x 4 C 2 



(32) 



2n; 



C - KC DX 4 C 2 - N B 



C -KC DX *C* -N B ) +Anl 



(33) 



At maximum value of n rnax — 3.694xl0 8 /xm -3 , the 
concentration of the substitutionally dissolved arsenic 
atoms C max is equal to 7. 943x10 s /im" 3 . The calculated 
curve n = n{C T ) agrees well with the experimental data, 
although only one fitting parameter KCd = 0.67xl0~ 15 
/jm 3 has been used. Full agreement is reached in the re- 
gion of the transition to saturation of the electron density 
n. This means that two arsenic atoms are incorporated in 
the cluster, and this arsenic cluster is doubly negatively 
charged, at least at the total dopant concentrations C T < 
1.0 x 10 9 ^m~ 3 . Minor differences in values of the theoret- 
ical curve and experimental data for C T > 2.3x 10 9 /j,m~ 3 
can arise due to heavy doping effects (change of the sil- 
icon zone structure, changes in the constants of forward 
Kl and backward Kr reactions, etc). It is clear from 
expressions (|22p. (j2"5)l that variation of the defect concen- 
tration Cd in the region of high impurity concentrations 
also influences the concentration of impurity atoms in- 
corporated in the clusters and, hence, the concentration 
of charge carriers. A change of Cd due to generation (ab- 
sorption) of the defects in the region with high arsenic 
concentration is quite possible. For example, such a gen- 
eration can occur due to the lattice expansion. Besides, 
it is possible that a high doping level leads to error in 
the determination of the experimental values for electron 
density, and the data for C T > 2.3xl0 9 /im~ 3 may be 
incorrect. 

Figs. 4 and b show the functions C = C(C* T ), C A = 
C A (C T ) 7 and n = n(C T ) calculated for clusters (DAS3) - 
and (DAs4)~ , respectively. As can be seen from these fig- 
ures, fitting to the experimental data is worse in compar- 
ison with Fig. 3b, especially for the total concentrations 
greater than 2.0 xlO 9 /im~ 3 . 

Note that multiple clustering can occur, too, and dif- 
ferent clusters can coexist in the equilibrium state. For 
example, the clusters incorporating different numbers of 
arsenic atoms (DAsi) - and (DAs 2 )~ or clusters in dif- 
ferent charge states (DAs 2 ) x , (DAs 2 ) _ , and (DAs 2 ) 2 ~ 
can coexist. As follows from the calculated curves for 
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FIG. 4: Calculated concentrations of substitutionally dis- 
solved arsenic atoms (solid line), clustered arsenic atoms 
(dashed line), and electron density (dotted line) against the 
total dopant concentration for the formation of singly nega- 
tively charged clusters incorporating three (a) and four (b) 
arsenic atoms. The measured electron density (circles) are 
taken from Solmi and Nobili* for the diffusion at a tempera- 
ture of 1050 °C. 



the electron density, when clusters (VAsi) - (Fig. 2) and 
(VAs 2 ) x (Fig. 1) are formed, it may be expected that the 
fall of the electron density in the plateau region will be 
less pronounced in the case of multiple clustering. How- 
ever, this problem requires further investigation. 



V. CONCLUSIONS 

Based on analysis of the experimental results and the- 
oretical calculations for clustering of arsenic atoms in sil- 
icon, it was shown that the developed clustering mod- 
els fail to explain the available experimental data, es- 
pecially the effect of electron density saturation at high 
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arsenic concentrations. Therefore, a new model of ar- 
senic clustering was proposed and analyzed. The main 
feature of the proposed model is the assumption that 
negatively charged arsenic complexes play a dominant 
role in the clustering process. To confirm this assump- 
tion, the concentration of the impurity atoms incorpo- 
rated in clusters and electron density were calculated as 
a function of the total arsenic concentration at a tem- 
perature of 1050 °C. Different cases of the formation of 
negatively charged clusters incorporating a point defect 
and one arsenic atom (DAsi) - , (DAsi) 2 ~ or more ar- 
senic atoms (DAs 2 )~, (DAs 2 ) 2 ~, (DAs 3 )~, (DAs 3 ) 2 ~, 
(DAs4)~, and (DAs4) 2 ~ were investigated. It was shown 
that in the case of (DAsi) - the concentration of charge 
carriers reaches saturation with increase in the total ar- 
senic concentration and concentration of the substitu- 
tionally dissolved impurity atoms. However, saturation 
was observed for very high values of the total arsenic 
concentration, conflicting with the experimental data. In 
the cases of doubly negatively charged cluster (DAsi) 2- 
and clusters incorporating more than one arsenic atom, 
the electron density reached its maximum value, then de- 



creased monotonically and slowly with an increase in the 
total arsenic concentration. The electron density calcu- 
lated for the formation of (DAs2) 2 ~ clusters agrees well 
with the experimental data and confirms the conclusion 
that two arsenic atoms participate in the cluster forma- 
tion. Minor difference between the theoretical curve and 
experimental values for C T > 2.3xl0 9 ^m -3 may be ex- 
plained by various heavy doping effects or uncertainty 
due to the experimental errors. 

Thus, the plateau on the experimental distribution of 
the charge carriers, characterized by a practically con- 
stant value of the electron density in the region of high 
doping level, may be attributed to the negatively charged 
clusters (DAs2) 2 ~. The parameters describing arsenic 
clustering at a temperature of 1050 °C were determined 
by fitting the calculated values of the electron density to 
the experimental data. The proposed model of cluster- 
ing by the formation of the negatively charged (DAS2) 2- 
system gives the best fit to the experimental data among 
all the present-day models and may be used in simulation 
of high concentration arsenic diffusion. 
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